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ABSTRACT 


A  theoretical  analysis  is  presented  for  the  problem  of  an 
axially  loaded  expansion  joint  or  bellow  consisting  of  a  shell  of 
revolution,  with  deep  periodic  corrugations.  The  problem  is  described 
analytically  by  two  simultaneous  differential  equations.  The  solu¬ 
tions  are  given  in  terms  of  functions  which  are  evaluated  by  numerical 
integration  and  tabulated. 

Experimental  results  are  presented  for  the  case  of  axial 
loading  which  agree  reasonably  with  the  theoretical  work. 
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INTRODUCTION 


The  use  in  piping  systems  of  expansion  joints,  or  bellows, 
consisting  of  thin  elastic  shells  of  revolution  is  of  importance  to 
the  design  engineer.  With  a  relatively  easy  method  of  determining 
the  stresses  and  displacements,  he  is  able  to  design  economically  for 
safety  as  well  as  for  performance. 

A  paper  by  Clark  [2],  in  1950,  presented  methods  of  asymptotic 
integration  for  the  analysis  of  thin  elastic  shells  of  revolution  loaded 
symmetrically  with  respect  to  their  axis.  Expansion  joints  subject 
to  axial  loads  and  corrugated  pipes  subject  to  axial  or  internal 
pressure  are  used  as  examples  of  the  general  procedure  developed  by 
Cl  ark . 

In  1962,  Laupa  and  Weil  [6]  presented  an  elastic  analysis 
of  U-shaped  expansion  joints  under  axial  loads  and  internal  or  ex¬ 
ternal  pressure.  Using  the  energy  method  for  the  toroidal  section 
and  combining  the  theory  of  flat  plates,  they  were  able  to  investigate 
any  U-shaped  bellows  under  any  arbitrary  combination  of  axial  and 
pressure  loading. 

More  recently,  Clark  [3]  has  solved  the  expansion  problem 
again  by  asymptotic  integration,  but  has  put  more  emphasis  on  the 
analysis  of  the  errors  involved  in  the  approximations.  The  cross- 
sectional  profile  he  considered  has  the  distinct  advantage  in  that 
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it  produces  a  smooth  surface  with  continuous  derivatives  such  that  both 
principal  curvatures  are  zero  at  a  point  of  inflection  of  a  meridian. 

The  analysis  presented  here  follows  the  method  used  earlier 
by  Reissner  [1]  to  solve  the  toroidal  shell  of  circular  cross-section. 

The  cross-sectional  profile  derived  by  Clark  [3]  is  introduced  into 
the  governing  differential  equations  and  a  solution  has  been  found 
through  numerical  integration  that  can  be  compared  with  the  numerical 
example  presented  by  Clark  [3]  and  with  the  experimental  results 
carried  out  by  the  author. 

In  the  solution  of  the  differential  equations,  the  terms 
inner  and  outer  solutions  are  used.  The  former  refers  to  that  region 
of  the  solution  where  the  boundary  layer  analysis  exists;  that  is, 
there  are  relatively  rapid  changes  in  the  dependent  variables.  The 
outer  solution  refers  to  that  portion  of  the  solution  which  does  not 
exhibit  a  boundary  layer  phenomenon.  In  this  region  the  bending  stresses 
as  well  as  the  stresses  due  to  rotation  of  the  shell  meridian  are  small 
compared  with  the  primary  membrane  stresses. 
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CHAPTER  I 

DERIVATION  OF  THE  EQUATIONS 

The  derivation  of  the  governing  equations  for  thin  elastic 
shells  of  revolution  follows  that  of  Reissner  [5].  The  same  notation 
that  he  used  has  been  adopted  here  and  is  used  throughout  the  thesis. 

Geometrical  Relations 

The  undeformed  middle  surface  of  a  shell  of  revolution  can 
be  described  in  the  form; 

r  =  r(£)  z  =  z(C) 

together  with  the  polar  angle  0.  The  angle  of  the  tangents  to  the 
meridian  (<j>)  is  given  by 


tan<j)  = 


dz 

dr 


(i.i) 


From  equation  (1.1),  the  following  relationships  can  be  derived: 


r 1 

cosd)  =  — 
a 


sincj)  = 


a 


(1.2) 


where  a  is  given  by 


4 


a 


o  1/2 

+  U')2} 


and  the  primes  denote  differentiation  with  respect  to 

The  principal  radii  of  curvature  of  the  middle  surface  of 
the  undeformed  surface  of  a  shell  of  revolution  are  given  by 


_L=  41 

a 


(1.3) 


1  _  si n4> 

R0  -  r  ' 


Analysis  of  Strain 

Denoting  the  undeformed  middle  surface  of  the  shell  by  the 
subscript  0,  the  deformed  middle  surface  can  be  written  in  the  form 

r  =  rQ  +  u  z  =  zQ  +  w  (1.5) 

where  u  and  w  are  the  displacements  in  the  radial  and  axial  directions 
respectively.  The  angle  <J>  of  the  deformed  surface  can  be  written 

(j)  =  4>0  -  3 

where  3  is  a  rotation  between  the  deformed  and  undeformed  meridians. 

The  assumption  is  made  here  that  the  deformations  due  to 
transverse  shear  strain  are  negligible  as  compared  to  one,  and  the 


r 


FIGURE  l.l  UNDEFORMED  6  DEFORMED 
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transverse  normal  stress  is  negligible  compared  with  the  remaining 
stresses . 

With  the  aid  of  figure  1.1  the  strains  along  the  centre 
line  of  the  shell  can  be  shown  to  be 


u 

£0M  rQ 

=  n  + 
coscj)  ^  r  1  ^ 

r  o 


(1.6) 

(1.7) 


where  the  subscript  M  refers  to  the  middle  surface. 

Equations  (1.3)  and  (1.4)  can  be  written  for  both  the  de¬ 
formed  and  the  undeformed  shell,  and  then  the  curvature  differences 
k  and  k0  are  given  by 


k,  =  -  ( — s_)  = 

5  “o  % 


Ke  =  '  (' 


sincj)  -  sin<f 


2-) 


(1.8) 


The  axial  displacement  w  can  be  found  by  substituting  equation 
(1.5)  into  equation  (1.2).  This  gives  the  axial  displacement  as 


(z  +  w)  =  asin<j) 

(1.9) 

w'  =  asin<t>  -  z  1 


The  compatibility  equation  is  derived  by  substituting  the 
radial  displacement  u  from  equation  (1.6)  into  equation  (1.7). 


,  !!!!o  n  +  (£en  . . 
e5M  cos<J>  r. ' 


which  reduces  to 


cos4>0  (rQ  e0fv])  -  cos<f>(r0'e  M)  =  r0'(cos<j)  -  cos<j>0)  (1.10) 

Equations  (1.7),  (1.8),  (1.9)  and  (1.10)  can  be  simplified  by 
expanding  in  powers  of  8,  and  neglecting  all  terms  except  the  first. 

The  radial  strain  now  becomes 


e 


y-r  -  etan<J)0 


fu '  tan<j> 


(1.11) 


the  circumferential  curvature  becomes 

8cos(J)o 


(1.12) 


the  axial  displacement  becomes 


w  =  z 


o 


-  (1  + 


(1.13) 


and  finally  the  compatibility  equation  becomes 
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cos ( roe0M)  '  (cos<t,0  +  6sin<i,o^ro'E5M'  =  ro'  (6sinV  U-14) 


Stress  Resultants  and  Couples 

Since  the  shell  is  axially  symmetric  and  only  axially  sym¬ 
metric  loads  are  being  considered,  there  will  be  no  dependence  upon 
the  polar  angle  0.  Because  of  the  independence  upon  0,  there  will  be 
only  four  stress  components  aQ,  a  ,  .  With  attention  being 

restricted  to  thin  shells,  terms  of  order  h/R  and  smaller  have  been 
neglected  as  compared  with  one,  and  o  is  neglected  compared  with  the 
other  stresses. 

Considering  the  normal  equilibrium  of  a  shell  element, 
(figure  1.2),  the  following  forms  of  the  stress  resultants  and  couples 
are  defined.  The  stress  resultants  are 


N 


5 


h/2 

» 

a  dc 

-h/2 


h/2 

N0  =  |  aQ  dC 
-h/2 


Q  = 


h/2 

/ 

-h/2 


d? 


and  the  couples  are 
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FIGURE  1.2  EQUILIBRIUM  OF  SMELL 
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h/2 

=  |  C  dc 

-h/2 

h/2 

M0  |  aQ  c  dc 

-h/2 


(1.15) 


For  convenience  the  same  notation  as  Reissner's  [4]  will  be 
employed  where  a  radial  (horizontal)  stress  resultant  H  and  an  axial 
(vertical)  stress  resultant  V  have  been  introduced.  These  are  related 
to  and  Q  in  the  following  way: 

=  Hcoscf>  +  Vs  i  n4>  (1.16a) 

Q  =  -  Hsin<f>  +  V cos4  .  (1.16b) 

Differential  Equations  of  Eguil ibriurn 

We  consider  the  equilibrium  in  the  axial  direction  as  shown 

« 

in  figure  1.2; 


(r(N^sincj)  +  Qcos<f>))  +  py  ra  =  0 

upon  substituting  equations  (1 . 16a)  and  (1 . 16b) 9  the  equilibrium  in  the 
axial  direction  is  given  by: 

I 

(rV)  +  p  ra  =  0  . 


(1.17) 


k 


\ 


. 


n 

Similarly,  the  equation  for  the  equilibrium  in  the  radial  direction 
and  the  equation  due  to  the  moments  about  a  ci rcumferenti al  line  are 


I 

(rH)  -  NgCt  +  p^|  ra  =  0 


(1.18) 


and 


(M^r)  +  ra(Hsin(J)  -  Vcos<}>)  -  M0acoscf>  =  0  (1.19) 

Stress-Strain  Relations 

Since  a is  neglected  as  compared  with  the  remaining  stresses, 
the  following  relationships,  derived  from  Hooke's  law, hold. 


where  E  is  Young's  modulus  of  elasticity  and  v  is  Poisson's  ratio. 
Inverting  these  equations,  the  stresses  can  be  written  in  terms  of 
the  strains  as  follows, 


(1.20) 
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The  strains  anywhere  in  the  shell  may  be  expressed  in  terms 
of  the  strain  along  the  middle  surface  plus  the  strain  due  to  the 
change  in  the  curvature.  The  strains  at  level  r,  are  given  by 


£0  =  £0M  +  CK0 


(1.21) 


Using  equations  (1.15),  (1.20)  and  (1.21),  it  can  be  easily 

shown  that 


Ce?M  -  vN0 

Ce0M  ~  N0  "  vN£ 
=  D(k^  +  vk0) 


(1.22) 


M0  =  D(k0  +  vk^) 

where 

Fh3 

C  =  Eh  D  =  - . 

12(l-v2) 

Equations  (1.22),  (1.19),  (1.18)  and  (1.17)  represent  a 
complete  system  of  seven  equations  that  can  be  solved  for  the  seven 


quantities  u,  3>  H,  V,  N0,  NU,  M0. 
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Equations  for  Small  Deflection 

The  small  deflection  theory  follows  from  the  differential 
equations  of  equilibrium,  together  with  the  equations  of  the  stress  re¬ 
sultants  and  by  neglecting  all  non-linear  terms  in  the  expressions  for 
the  strains.  The  resulting  system  of  equations  can  be  reduced  to  two 
simultaneous  second  order  differential  equations  with  the  basic  vari¬ 
ables  3  and  (r  H).  In  terms  of  these  variables,  the  system  of  equations 
become , 


r  V  =  -  r  a  p  d£ 

0  0  0  V 

* 

(1.23) 

rQ  =  (rQ  H)  cos  <J>Q  +  (rQ  V)  sin  4>0 

(1.24) 

rQ  Q  =  -  (rQ  H)  sin  +  (rQ  V)  cos  <J>Q 

0-25) 

a  N„  =  (r  H)'  +  r  a  pu 

00  v  o  o  o 

(1.26) 

aQ  =  D ( 3 '  +  v(ro7rQ)3) 

(1.27) 

aQ  M0  =  D((r0'/r0)3  +  v3') 

(1.28) 

u  =  (r0/c)(Ne-vN^ 

(1.29) 

w  =  |  {(Zq'/C) (N^-vNq)  -  rQ ' / 3 }  d£ 

(1.30) 

The  first  of  the  two  simultaneous  differential  equations  for 
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3  and  r  H  is  obtained  by  introducing  the  moments  Mr  and  M.  from  equations 
(127)  and  (l  .28)  into  the  moment  equilibrium  equation  (1 . 19)  for  the  middle 
surface.  Assuming  that  the  stiffness  factors  C  and  D  are  constant  the 
equation  becomes. 


3 


+ 


(r  /a  ) ' 
o _ oy 

(r  /a  T~ 

v  o  o' 


3'  -  {(zr-) 


-  v 


(ro,/ao)']  6 


z  1  (r  H)  r  '  (r  V) 

( r  D/a  )  (r  D/a  ) 

v  o  o  o  o 


In  a  similar  manner,  the  second  differential  equation  is  de¬ 
rived  by  introducing  the  strains  and  from  the  set  of  equations 
(1.22)  along  with  and  NQ  from  equations(  1  .24)and  (l  .26)  into  the  com¬ 
patibility  equation  (1.14). 


(  K*  /  (X  )  * 

(r0  H)"  +  Tk/%T  (ro  H)’  -  {(ro7ro) 


(r  7a  )'  z 

{rjaj  o  •  (r0/a0  c) 


o  o' 


{[rr°7a°r  *  v(ro'/ro)]  (ro  “o  pl4 


+  (ro  %  PH),}  +  {V  fz  +  Tr°7aJ'"'}  (ro  V) 

r  0 


0  O' 


+ 


v(zo'/r0)  (r0  V) 


II 


Equations  I  and  II  completely  describe  the  behavior  of  any 


15 

symmetric  shell  of  revolution  where  the  theory  of  small  deflections  is 
val id. 

The  type  of  cross-section  to  be  considered  is  shown  in  figure 
1.3.  Clark  [3]  has  used  the  parametric  equations 

r  =  a  +  b  sin  ^ 

and  z  =  -  c  (C  -  sinC) 

where  c"  =  I  Cl  <  S- 

7T-  2  11  —  2 

to  describe  the  profile  between  the  points  A  and  B. 

Figure  1.4  shows  a  comparison  between  the  cross  section  used 
here  and  that  of  a  U  shape  profile  for  the  case  b  =  2c.  The  difference 
between  the  two  profiles  seems  to  be  negligible  as  compared  to  the 
variations  in  the  shape  of  a  manufactured  bellows. 

When  the  only  load  on  the  shell  is  an  axial  edge  load,  the 
overall  equation  of  axial  equilibrium  shows  that  V  satisfies  the  equation 

2tt  r  V  =  P 
u 

where  P  is  the  total  axial  load  and  since  there  is  assumed  to  be  no 
internal  or  external  pressure  py  =  p^  =  0. 

Substituting  the  quantities  r  ,  zQ  and  aQ  into  the  governing 
differential  equations  I  and  II,  it  can  be  verified  that  they  reduce  to 
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FIGURE  1.3  CROSS  -  SECTIONAL  PROFILE 
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FIGURE  1.4  ASSUMED  SHAPE  COMPARED 

TO  A  CIRCULAR  ARC 
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the  forms. 


2  2 
oi.  ,  f  Xcos^  sing(cos^-n  (1-cosQ)-, p,  f(  Acos^  % 

B  1+Asi"?  "(cosW(l-cos5)2)  '  1+Asin« 


,  r  Xsinr  Acosgsing(cosg-ri  (l-cosg)) 
"1+Asin«  "  (cos2  £+n^(  l-cosq)^)(l+Xsin^) 


cb  (1-cosg) (cos^+r^( 1-cosg)^) 
aD  ~  (1+Asini;) 


1/2 


(roH) 


cos£;(cos^+ri^(  1  -cosg)^) 
aD  /l+Asin^T” 


1/2 


( roV) 


(1.31) 


and 


(rH)"  +  ^TWsfnFT  -  gMLcos|t.nJl-c°sO.).}  ( r  H) 
o  ( 1+Asinf,;  (cos25+o2(1-cos5)  )  0 


2  2 
r /  Acosg  v  vAsin^  vAcos^si n^(-cosg+n  (1-cosg))  /  „ \ 

' l+Asi n?  -  l+Asin?  -  (1+As1nC)( 

cos^+n2  (1-COS5)2)  0 


2  2  2  1/2  - 

+  bcC  (  1  -COSg)  (cOS^-tT]  (  1  -COSg)  J  _  £  ,--A(  1  -cosQcosg; 


( l+Asin^) 


(1+Asin^)1 


+  v  {  -- 


sin^  |  y  1-LU^p  I  ii^y-uuac,T[|  V  I  -mac,;  j  j  j  (  r  V)  (1  32) 

fl+Asiniy  (1+AsinC) (cos2^+n2  (l-cos£)2)  0 


(1-cos£;)sing(-cosC+Ti  (1-cosg)) 


where  A  and  n  have  been  defined  as 


Equations  (1.31)  and  (1.32)  can  be  solved  subject  to  appropriate 


boundary  conditions.  The  four  boundary  conditions  are  due  to  the  sym¬ 
metry  conditions  of  vanishing  3  and  Q  for  £  =  tt/2.  These  conditions 
can  be  written  in  the  form; 


20 


CHAPTER  II 

SOLUTION  OF  EQUATION 


2 . 1  Outer  Sol ution 

To  determine  to  what  extent  the  solution  of  equations  (1 . 31) 
and  (l  .32) does  not  exhibit  a  boundary  layer  phenomenon,  the  dependent 
variables  are  non-dimensional i zed  in  the  following  way, 


B  =  B0  FU)  (r0H)  =  >PoG(5)  (rQV)  =  ^  (2.1) 


where 


(2.2) 


Using  these  non-dimensional i zed  variables,  the  two  governing 
di  fferential  equations  (1 .31)  and  (1  .32)  become. 
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Since  we  are  free  to  choose  3  ,  the  coefficients  of  the 


coupling  terms  are  equated  to  a  new  parameter  y 
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From  (2.2)  and  (2.5),  it  follows  that 
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The  differential  equations  (2.3)  and  (2.4)  now  become 
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(2.9) 


with  the  boundary  conditions 


at  E,  -  +_  ^  ,  F  -  0  G  -  0 


Expressions  for  the  stress  resultants  and  stress  couples,  in 
terms  of  F(£)  and  G(£)  are  obtained  from  the  solution  of  equations 
(2.8)  and  (2.9) 


N5  =  {Gcos5+ sin5}  (2-10) 


N.  = - o  ^2 - H/2  G' 

2irb(cos£5+n  (1-cosC)  ) 


(2.11) 


23 


NL  = 


PnX 


2Try(cos2£+r|2(l--cos02) 


(F1 


,  vXcos^ 

+  THxil'r ifT  F> 


(2.12) 


M 


0 


Pnx 

27iy(  cos2^+n2  (l-cosC)2)1^2 


Acosg 
1+Asi n£ 


F  +  vF'}  . 


(2.13) 


For  a  solution  to  equations  (2.8)  and  (2.9)  to  exist,  valid 
throughout  the  interval  -  tt/2  <  £  <  tt/2,  the  parameter  y  would  have  to 
be  of  order  of  unity.  This  can  be  seen  from  equation  (2.8),  for  as  y 
becomes  large  compared  to  unity,  the  asymptotic  solution  for  G(£)  is 
given  by 


GU)  ~ 


COSE, 

1 -cos^ 


which  has  a  singular  point  at  £  =  0.  With  y  of  order  of  unity,  A  must 
be  small  as  compared  to  unity  since  from  equation (2 . 7)  A/y  =  ( h/c ) / [1 2 ( 1  -v  )]  , 

which  is  always  small  when  compared  to  unity.  Neglecting  terms  in  A 
as  compared  to  the  remaining  terms,  equations  (2.8)  and  (2.9)  become 
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cos^+ri  ( 1 -cos£; 

1/2 

+  y(l-cos^)  (cos2^+n?'(l-cos^)2)  F  =  0  (2.15) 

together  with  the  boundary  conditions 

at  ?  =  F  =  0  G  -  0 

It  should  be  noted  here,  that  as  y  becomes  large  compared 
with  unity,  the  equations  are  no  longer  valid  near  the  origin.  In 
addition,  although  A/y  is  always  small  compared  with  one,  the  value 
of  A  is  no  longer  necessarily  small. 

For  an  approximation  to  the  outer  solution  for  the  case 

b/a  =  0.2117  c/b  =  0.6257  h/a  =  0.0069 

i 

equations  (2.14)  and  (2.15)  have  been  solved  numerically  with  F(£)  and 
G(C)  plotted  in  figures  2.1  and  2.2.  This  example  corresponds  to  the 
experimental  bellow  that  is  used  in  Chapter  III.  The  geometric  values 
given  above  correspond  to  the  case  y  =  11.7640,  A  =  0.2117. 

The  radial  stresses  are  shown  in  figure  2.3  and  the  circum¬ 
ferential  stresses  in  figure  2.4  for  a  compressive  load.  The  values 
are  not  valid  near  6  =  0>  but  become  closer  to  the  actual  stresses  as 
£  -*  +  If.  The  outer  solution  obtained  here  will  be  matched  up  with 
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the  inner  solution  to  give  a  more  accurate  solution  throughout  the 
range  of  ?. 

2.2  Inner  Solution 

From  the  structure  of  the  differential  equations  (2.8)  and 
(2.9),  it  can  be  seen  that  when  p  is  large  compared  to  unity,  the 
neighborhood  of  ?  =  0  is  an  area  of  rapid  changes  in  the  dependent 
variables  F(?)  and  G(?).  For  the  analysis  of  the  inner  solution  a 
new  independent  variable  p  is  introduced  in  the  form 

P  =  T? 

where  x  »  1  such  that  a  narrow  region  surrounding  ?  =  0  becomes  a 
region  surrounding  p  =  0.  Since  we  are  considering  the  neighborhood 
of  ?  =  0  the  approximations 

2 

sin?  =  —  cos?  =  1  - 

T  2x2 


will  be  used. 

Defining 


F  =  Fq  f(p)  G  =  Gq  g(p) 


and  using  the  above  approximations,  the  governing  differential  equations 
(2.8)  and  (2.9)  become 
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Two  conditions  for  the  determination  of  the  parameters  t,  Fq 
and  Gq  come  from  equation  (2.16),  by  requiring  that  the  coefficients  of 
f ' '  ,  g  and  1  be  of  comparable  magnitude 
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The  third  condition  follows  from  equation  (2.17)  by  requiring  that 
the  terms  with  g1'  and  f  be  of  comparable  magnitude 


(2.19) 


Solving  equations  (2.18)  and  (2.19)  for  the  parameters  t, 


F  and  G  ,  we  have 
o  o 
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Using  the  relationships  (2.20)  to  (2.22)  and  the  fact  that  t  »  1,  the 
two  simultaneous  differential  equations  (2.16)  and  (2.17)  can  be  reduced 
by  neglecting  terms  of  smaller  magnitudes 


-  p2g  =  1 

(2.23) 

+  P2f  =  0 

(2.24) 

with  the  boundary  conditions  in  the  form 


at  p  =  ±  w  ,  f=0  g  =  0 


It  must  be  noted  that  the  solution  of  equations  (2.23)  and  (2.24)  are 
val id  only  if  t  »  1 . 

Expressions  for  the  stress  resultants  and  couples  follow 
from  equation  (2.10)  to  (2.13)  in  the  form, 
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The  total  stress  in  the  surface  fibers  of  the  shell  can  be 
written  as  the  sum  of  the  direct  membrane  stress  and  the  stress  due  to 
bending,  o  =  .  The  upper  sign  refers  to  the  surface  which  is 

farthest  from  the  axis  of  revolution,  while  the  lower  sign  refers  to 
the  surface  closest  to  the  axis  of  revolution.  The  direct  stresses 
are  related  to  the  stress  resultants  by 

N  N 

a£D  ~  'IT  a0D  =  T  ('2,29 

while  the  bending  stresses  are  related  to  the  stress  couples  by 


(2.30) 


Using  equations  (2.25)  to  (2.30),  the  direct  stress  and  bending  stress 
for  both  the  radial  and  circumferential  directions  can  be  written  in 
the  form, 
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Equations  (2.23)  and  (2.24)  are  solved  by  numerical  integration  and 
the  functions  f(p)  and  g(p)  are  plotted  in  figures  2.5  and  2.6  re¬ 
spectively.  Since  the  solution  of  the  equations  is  independent  of  y, 
other  than  y  being  large  as  compared  to  unity,  the  values  of  f(p)  and 
g(p)  are  tabulated  in  appendix  C  to  be  used  in  determining  the  stresses 
in  the  region  near  £  =  0  for  all  U-shaped  bellows. 


Using  the  example  cited  earlier  in  the  outer  solution,  the 


stress  distribution  is  solved  for  the  inner  solution.  The  radial 
direct  stress,  the  radial  bending  stress  and  the  total  radial  stress 
for  the  inner  surface  are  plotted  in  figure  2.7.  Similarly  in  figure 
2.8,  the  ci rcumferenti al  direct  stress,  the  circumferential  bending 
stress  and  the  total  circumferential  stress  are  plotted. 
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CHAPTER  III 

EXPERIMENTAL  INVESTIGATION 

A  manufactured  Flexonics  expansion  joint  was  purchased  for 
the  experimental  investigation  of  the  stresses  caused  by  an  axial  load. 
The  bellow  pictured  in  figure  3.1  is  designed  to  fit  a  nominal  pipe 
size  of  eight  inches.  With  three  corrugations ,  it  has  a  maximum  axial 
traverse  of  one  and  a  half  inches.  The  bellow  is  manufactured  out  of 
stainless  steel  type  304  and  is  free  flexing  in  that  there  is  no  control 
on  the  maximum  or  minimum  traverse  travel  or  angular  rotation.  The 
actual  contour  of  the  corrugations  as  shown  in  figure  3.3,  is  more 
of  an  S-shape  profile  rather  than  the  shape  used  earlier  in  the  theo¬ 
retical  analysis. 

To  determine  the  values  of  a,  b,  c,  and  h  for  the  bellow  con¬ 
sidered  here  a  number  of  readings  were  taken  and  the  average  values 
used.  These  readings  are  recorded  in  appendix  A  where  the  average 
values  are  given.  A  variation  of  0.065  inches  was  found  in  the  depth 
while  a  variation  of  0.040  inches  was  found  in  the  width  of  the  corru¬ 
gations.  The  readings  for  the  thickness  of  the  material  were  taken 
on  the  end  flanges  using  a  micrometer. 

Four  Budd  metal  film  strain  gauges  (type  C6-121)  were  placed 
at  ninety  degrees  to  each  other  on  each  of  the  three  corrugations.  The 
reason  for  this  arrangement  was  to  try  to  predict  if  a  purely  axial 
load  was  being  applied.  Since  the  loads  being  applied  were  in  the 
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range  from  zero  to  five  hundred  pounds,  the  strain  at  the  tips  of  each 
corrugation  was  too  small  to  give  any  indication  as  to  whether  an 
eccentric  load  was  being  applied.  To  determine  the  stress  distribu¬ 
tion  along  a  meridian,  fourteen  strain  gauges  were  mounted  on  the  out¬ 
side  surface  of  one  half  of  the  middle  corrugation  as  shown  in  figure 
3.2.  Grouping  the  gauges  into  pairs,  one  to  record  the  radial  strain 
and  one  to  record  the  circumferential  strain,  and  equally  spacing  each 
pair  along  a  meridian,  the  complete  strain  distribution  was  determined. 
Opposite  the  pairs  of  gauges  on  the  outside  surface  were  mounted  an¬ 
other  five  pairs  of  gauges  to  record  the  strain  distribution  on  the 
inside  surface.  In  figure  3.4,  the  location  of  all  the  strain  gauges 
is  shown  and  the  numbers  correspond  to  the  gauges  in  appendix  B. 

An  Instron  testing  machine  was  used  for  the  application  of 

s 

the  axial  load  as  shown  in  figure  3.5.  The  maximum  load  range  used 
was  five  hundred  pounds  with  an  accuracy  of  plus  or  minus  two  pounds. 
Load  increments  of  one  hundred  pounds  were  applied  and  the  strains 
were  recorded  by  using  two  types  of  strain  indicators  -’  a  Baldwin- 
Lima-Hami 1  ton  strain  indicator  and  an  Automation  Industries  strain 
indicator.  Gauges  one  to  seven  were  hooked  to  the  former  while  gauges 
eight  to  twenty-four,  using  two  switching  units,  were  hooked  to  the 
latter.  All  experimental  data  has  been  listed  in  appendix  B. 

From  equation  (1.20),  the  radial  and  circumferential  stresses 
can  be  determined.  For  comparison  purposes,  the  stresses  are  written 
in  the  dimensionless  form 
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FIG.  3.5  TEST  SET-UP 
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The  values  ch  o^/P  for  the  outside  radial  gauges  are  listed 
in  table  3.1.  Similarly  for  the  outside  circumferential,  the  inside 
radial,  and  the  inside  circumferential  gauges,  the  values  of  cha/P  are 
listed  in  tables  3.2,  3.3  and  3.4  respecti vely . 


45 


Table  3.1 
Outside  Radial 


£(rads) 

ch 

a^/P 

Average 

tt/2 

0.25 

0.25 

0.25 

0.25 

0.25 

0.9407 

1 .28 

1.45 

1.50 

1.50 

1.43 

0.4555 

1.97 

2.06 

2.04 

1.97 

2.01 

0 

-0.12 

-0.04 

-0.05 

-0.10 

-0.08 

-0.4555 

-0.87 

-1.09 

i 

_ t 
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VO 

C\J 
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• 

1 

-1.02 

-0.9407 

-2.34 

-2.53 

-2.47 

-2.42 

-2.44 

-tt/2 

-0.74 

-0.88 

-0.93 

-0.97 

-0.88 

Table  3.2 

Outside  Circumferential 

C(rads) 

ch  a0/P 

Average 

tt/2 

0.83 

0.85 

0.85 

0.84 

0.84 

0.9407 

0.83 

0.87 

0.86 

0.84 

0.85 

0.4555 

-0.28 

-0.29 

-0.31 

-0.34 

-0.31 

0 

-1  .18 

-1.16 

-1.13 

-1.13 

-1  .15 

-0.4555 

-0.14 

-0.17 

-0.15 

-0.12 

-0.15 

-0.9407 

-0.40 

-0.43 

-0.39 

-0.37 

-0.40 

-tt/2 

-0.93 

-0.95 

-0.94 

-0.94 

-0.94 

t 
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Table  3.3 
Inside  Radial 


^(rads) 

ch  o £. 

/P 

Average 

-tt/2 

0.90 

0.94 

0.96 

0.96 

0.94 

-0.9407 

2.34 

2.34 

2.27 

2.21 

2.29 

-0.4555 

1.51 

1.50 

1.44 

1.39 

1.46 

0 

-0.16 

-0.06 

0.02 

0.08 

-0.03 

0.4555 

-1.74 

-1.65 

LO 

r— 

1 

-1.44 

-1.59 

Table  3.4 

Inside  Circumferential 

C(rads) 

ch  aQ/P 

Average 

-tt/2 

-0.42 

-0.39 

-0.36 

-0.34 

-0.38 

-0.9407 

1.08 

1.11 

1.10 

1.10 

1.10 

-0.4555 

1 .38 

1 .38 

1  .34 

1.31 

1.35 

0 

-0.54 

-0.54 

-0.52 

-0.50 

-0.53 

0.4555 

-1.11 

-1.10 

-1.08 

-1.06 

-1.09 
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CHAPTER  IV 

DISCUSSION  OF  RESULTS 

The  complete  theoretical  stress  distribution  is  found  by 
matching  the  inner  and  outer  solutions.  Figures  4.1  and  4.2  show 
the  two  solutions  superimposed  for  the  radial  and  circumferential 
stresses  respecti vely .  The  complete  solution  can  be  broken  down  into 
two  regions,  one  where  only  the  inner  solution  is  valid,  and  one  where 
only  the  outer  solution  is  valid. 

In  the  interval  £  -  ±  0.40,  the  deviations  of  the  outer 
solution  from  the  inner  solution  become  negligible.  This  is  more 
likely  a  coincidence  since  for  large  values  of  y,  the  outer  solution 
is  not  valid  near  5=0.  As  5  increases  from  +0.40  to  +_  1 .20,  the 
variation  between  the  two  solutions  increases  to  a  maximum  near 
5  =  +0.65. 

Throughout  this  region  the  inner  solution  is  valid,  while 
for  5  greater  than  +_  1.20,  the  outer  solution  becomes  valid.  The 
matching  of  the  inner  and  outer  solutions  produces  a  relatively  con¬ 
tinuous  solution  for  the  stress  distribution  in  the  expansion  joint. 

For  a  qualitative  view  of  the  theory  developed  here,  the 
experimental  values  from  Chapter  III  and  the  results  of  Clark  [3] 
are  used  for  comparison. 

In  figure  4.3,  both  the  inside  and  the  outside  radial  stress 
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FIGURE  4.1  INNER  S  OUTER 

SOLUTIONS  (RADIAL) 
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FIGURE  4.2  INNER  8s  OUTER 

SOLUTIONS  (CIRCUMFERENTIAL) 
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distributions  are  plotted  to  show  the  difference  between  the  theoretical 
and  experimental  results.  Since  the  direct  radial  stress  is  small 
as  compared  to  the  bending  stress,  the  difference  in  magnitude  between 
the  inside  and  outside  stress  is  negligible.  In  the  region  surrounding 
the  origin,  the  present  solution  follows  that  of  Clark's  very  closely. 
However,  for  E,  greater  than  +_0.40  the  variation  becomes  very  distinct. 
At  £  =  0.85,  the  present  solution  predicts  a  stress  of  25%  higher  than 
that  of  Clark's.  For  larger  values  of  £,  the  two  inner  solutions  con¬ 
verge  to  the  outer  solution.  The  solution  presented  here  is  in  better 
agreement  with  the  experimental  points  than  Clark's,  thus  indicating 
the  validity  of  the  solution  in  the  region  of  maximum  stress. 

The  inside  and  outside  experimental  stresses  and  the  cir- 
cumferential  stress  for  both  theories  are  plotted  in  figure  4.4. 

The  difference  between  the  theories  is  more  pronounced  throughout 
the  range  of  £  than  was  the  case  for  the  radial  stress.  Again  at  the 
point  of  maximum  stress,  the  present  theory  predicts  a  circumferential 
stress  of  25%  higher  than  that  by  Clark.  Both  theories  converge  to 
the  outer  solution  for  large  values  of 

The  stress  distributions  in  figures  4.3  and  4.4  are  plotted 
for  the  inside  surface  of  a  bellow  under  a  compressive  axial  load. 

The  stresses  on  the  inside  surface  give  a  stress  distribution  similar 
to  the  theory.  For  the  outside  circumferential  stresses,  the  stress 
corresponding  to  N0  is  of  opposite  sign  to  the  bending  stress,  yielding 
a  distribution  completely  different  from  that  of  the  inside  surface. 
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FIGURE  4.4  CIRCUMFERENTIAL 

STRESSES  INSIDE  SURFACE 
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The  poor  agreement  between  the  experimental  and  the  theoretical  values 
for  the  circumferential  stress  on  the  outer  surface  can  be  contributed 
to  the  variation  of  the  wall  thickness  throughout  a  corrugation. 

Another  source  of  error  would  be  due  to  the  difficulty  in  mounting 
the  strain  gauges;  however,  the  present  theoretical  analysis  seems  to 
fit  reasonably  well  to  the  experimental  results. 
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CHAPTER  V 

CONCLUDING  REMARKS 

The  major  conclusions  that  can  be  drawn  from  the  theoretical 
analysis  and  experimental  investigation  are: 

(1)  The  stresses  are  very  sensitive  to  variations  in  the  wall 
thickness.  An  increase  of  0.005  inches  in  h  produces  an  increase  in 
the  experimental  stress  by  15%,  while  the  theoretical  stress  decreases 
by  15%. 

(2)  Because  of  the  high  dependence  upon  h,  manufacturers  should 
incorporate  in  their  specifications  the  minimum  wall  thickness  of 
the  bellow. 

(3)  The  experimental  and  theoretical  values  calculated  for  the 
stresses  do  not  differ  by  more  than  30%.  Since  an  average  value  of 
.035  inches  was  used  for  the  bellow  thickness,  and  the  actual  values 
vary  from  this  figure  up  to  .005  inches,  the  discrepancy  between 
theoretical  and  experimental  stress  values  is  within  the  bounds 

of  experimental  error. 
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APPENDIX  A 


DETERMINATION  OF  EXPERIMENTAL  PARAMETERS 

Average  values  for  the  physical  parameters  a,  b,  c  and  h 
were  determined  prior  to  any  experimental  tests.  Vernier  calipers 
were  used  for  all  measurements. 


h  =  material  thickness 
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Measurements 


Average 


Depth 
(i nches) 

2.157 
2.182 
2.165 
2.130 
2.177 
2.168 
2.165 
2.190 
2.192 
2.160 
2.127 

2.158 


Wi  dth 
( i nches) 

1.365 

1.363 

1.330 

1.360 

1 .325 

1  .353 

1.342 

1.335 

1  .330 

1.340 

1 .325 

1  .341 


Using  the  average  values  from  above,  the  particular  bellow 
considered  here  has  the  following  dimensions 


a  -  5.074" 
b  =  1.074" 
c  =  0.672" 


h  =  0.035" 


APPENDIX  B 


EXPERIMENTAL  DATA 


Test  #1 


60 


Strain  Gauge 

(micro  inches  Load  (pounds) 


per  inch) 

100 

200 

300 

400 

1 

2 

135 

330 

525 

716 

3 

285 

605 

902 

1185 

4 

34 

102 

146 

170 

5 

-108 

-270 

-420 

-552 

6 

-290 

-638 

-960 

-1260 

7 

-  45 

-142 

-250 

-350 

8 

100 

210 

320 

434 

9 

60 

120 

172 

220 

10 

-110 

-242 

-375 

-510 

11 

-145 

-306 

-458 

-600 

12 

15 

42 

70 

102 

13 

40 

90 

145 

200 

14 

-  86 

-182 

-266 

-352 

15 

-  83 

-175 

-260 

-340 

16 

50 

no 

174 

244 

17 

122 

248 

375 

494 

18 

-  60 

-134 

-208 

-280 

19 

-  68 

-158 

-246 

-337 

20 

130 

282 

433 

582 

21 

256 

538 

794 

1036 

22 

136 

288 

426 

548 

23 

-  6 

18 

62 

123 

24 

-182 

-362 

-504 

-628 

Test  #2 
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Strain  Gauge 

(micro  inches  Load  (pounds) 

per  inch)  100  200  300  400 


1 


2 

166 

360 

571 

733 

3 

306 

630 

945 

1195 

4 

54 

106 

165 

155 

5 

-  95 

-275 

-402 

-572 

6 

-292 

-659 

-960 

-1290 

7 

-  93 

-150 

-235 

-380 

8 

109 

220 

327 

434 

9 

64 

124 

175 

222 

10 

-123 

-257 

-392 

-530 

11 

-163 

-325 

-477 

-624 

12 

16 

46 

76 

112 

13 

42 

92 

148 

204 

14 

-101 

-196 

-280 

-364 

15 

-100 

-194 

-277 

-357 

16 

54 

118 

180 

248 

17 

130 

264 

386 

510 

18 

-  70 

-150 

-222 

-295 

19 

-  85 

-170 

-263 

-355 

20 

148 

304 

455 

606 

21 

287 

572 

826 

1072 

22 

158 

312 

446 

570 

23 

5 

28 

76 

133 

24 

-201 

-375 

-516 

-637 

Test  #3 


62 


Strain  Gauge 

(micro-inches  Load  (pounds) 

per  inch)  100  200  300  400 


1 


2 

153 

330 

513 

706 

3 

296 

600 

894 

1170 

4 

33 

66 

96 

120 

5 

-124 

-315 

-480 

-606 

6 

-329 

-705 

-1038 

-1340 

7 

-  56 

-190 

-309 

-414 

8 

111 

222 

330 

436 

9 

64 

124 

174 

221 

10 

-130 

-263 

-400 

-540 

11 

-170 

-332 

-486 

-633 

12 

18 

46 

76 

108 

13 

43 

95 

148 

205 

14 

-106 

-200 

-288 

-372 

15 

-104 

-198 

-280 

-363 

16 

60 

121 

184 

250 

17 

139 

267 

392 

504 

18 

-  76 

-152 

-230 

-303 

19 

-  90 

-180 

-270 

-364 

20 

154 

310 

461 

614 

21 

300 

582 

840 

1084 

22 

164 

318 

450 

574 

23 

0 

32 

81 

135 

24 

-206 

-382 

-524 

-642 

Test  #4 
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Strain  Gauge 

(micro-inches  Load  (pounds) 

per  inch)  100  200  300  400 


1 


2 

131 

326 

500 

687 

3 

278 

594 

880 

1150 

4 

12 

65 

90 

101 

5 

-145 

-320 

-479 

-624 

6 

-347 

-713 

-1040 

-1355 

7 

-  70 

-190 

-306 

-429 

8 

110 

222 

330 

436 

9 

63 

124 

174 

221 

10 

-130 

-264 

-400 

-540 

11 

-170 

-333 

-485 

-633 

12 

18 

46 

75 

107 

13 

42 

93 

148 

205 

14 

-105 

-201 

-288 

-371 

15 

-106 

-200 

-284 

-362 

16 

53 

120 

181 

252 

17 

134 

269 

388 

510 

18 

-  75 

-154 

-228 

-303 

19 

-  90 

-177 

-271 

-366 

20 

152 

310 

460 

612 

21 

302 

584 

838 

1082 

22 

162 

318 

450 

572 

23 

2 

38 

79 

137 

24 

-206 

-382 

-524 

-644 

Average  of  Tests 
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Strain  Gauges 
(mi cro-i nches 
per  inch) 

100 

Load  (pounds) 

200  300 

400 

1 

2 

146 

336 

527 

710 

3 

291 

607 

905 

1175 

4 

33 

85 

124 

136 

5 

-118 

-295 

-445 

-558 

6 

-314 

-679 

-999 

-1311 

7 

-  66 

-168 

-276 

-393 

8 

107 

218 

327 

435 

9 

63 

123 

174 

221 

10 

-123 

-256 

-392 

-530 

11 

-162 

-324 

-476 

-622 

12 

17 

45 

74 

107 

13 

42 

92 

147 

203 

14 

-100 

-195 

-280 

-365 

15 

-  98 

-192 

-275 

-355 

16 

54 

117 

180 

248 

17 

131 

262 

385 

504 

18 

-  70 

-147 

-222 

-295 

19 

•  -  83 

-171 

-262 

-355 

20 

146 

301 

452 

603 

21 

286 

569 

824 

1068 

22 

155 

309 

443 

566 

23 

0 

29 

74 

132 

24 

-199 

-375 

-517 

-638 
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APPENDIX  C 

VALUES  OF  f(p)  AND  g(p) 

Calculations  of  the  function  f  and  g 

The  numerical  method  used  consisted  of  a  Hammings  modified 
predictor  corrector  for  the  integration  of  equations  (2.23)  and  (2.24). 
In  order  to  use  this  procedure,  the  initial  conditions  f(0)  ,  f'(0) 
g (0)  ,  g 1 ( 0 )  must  be  specified.  We  know  two  conditions  since  the 
functions  f  and  g  are  both  even  values  in  p.  The  other  initial  con¬ 
ditions  were  found  by  adapting  Newton's  method  for  quadratic  poly¬ 
nomial  factors  [7].  The  criterion  used  for  the  determination  of  the 
correct  value  was  that  the  functions  should  be  equal  to  their  asymptotic 
expansion  at  p  =  5.0.  Since  only  one  condition  can  be  satisfied  for 
each  function  it  is  impossible  to  specify  that  the  slopes  of  the 
functions  should  be  equal  to  those  of  their  asymptotic  expansions. 
However,  in  the  region  p  =  5.0,  the  slopes  do  match  within  the  accuracy 
specified  for  the  values  of  the  functions. 

In  figure  2.5  and  2.6,  the  values  of  f  and  g  are  found  by 
the  numerical  integration  for  p  <  5.0,  while  for  p  >  5.0,  the  asymptotic 
expansion  is  used.  These  values  of  f ,  f ' ,  f  and  g'  are  similar  to  the 
function  T^  g(x)  which  is  described  in  [3],  The  relations  between 
these  functions  are 
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f(p)  = 

-  Re 

T2,0 

(p) 

g(p)  = 

Im 

T2  ,0 

(P) 

f(p)  = 

-  Re 

— 1 

ro 

s* 

o 

(P) 

g1 (p)  = 

Im 

T  1 

‘2,0 

(p) 

The  values  of  the  functions  and  their  derivatives  for  p  < 
given  below. 


5.0  are 
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p 

f(p) 

ftp) 

g(p) 

g'(p) 

0 

-0.9270 

0.0000 

-0.9270 

0.0000 

0.1 

-0.9220 

0.0997 

-0.9270 

0.0003 

0.2 

-0.9072 

0.1975 

-0.9269 

0.0024 

0.3 

-0.8827 

0.2917 

-0.9264 

0.0081 

0.4 

-0.8490 

0.3802 

-0.9251 

0.0188 

0.5 

-0.8069 

0.4615 

-0.9225 

0.0356 

0.6 

-0.7570 

0.5336 

-0.9178 

0.0593 

0.7 

-0.7005 

0.5948 

-0.9104 

0.0901 

0.8 

-0.6385 

0.6439 

-0.8995 

0.1278 

0.9 

-0.5722 

0.6793 

-0.8846 

0.1715 

1  .0 

-0.5031 

0.7003 

-0.8651 

0.2199 

1.1 

-0.4326 

0.7062 

-0.8405 

0.2714 

1.2 

-0.3624 

0.6969 

-0.8108 

0.3239 

1.3 

-0.2937 

0.6730 

-0.7758 

0.3750 

1.4 

-0.2282 

0.6351 

-0.7359 

0.4224 

1.5 

-0.1671 

0.5850 

-0.6915 

0.4637 

1 .6 

-0.1116 

0.5247 

-0.6434 

0.4969 

1.7 

-0.0624 

0.4566 

-0.5925 

0.5203 

1  .8 

-0.0204 

0.3832 

-0.5397 

0.5328 

1.9 

0.0141 

0.3077 

-0.4863 

0 . 5336 

2.0 

0.0411 

0.2330 

-0.4334 

0.5226 

2.1 

0.0609 

0.1620 

-0.3821 

0.5009 

2.2 

0.0737 

0.0969 

-0.3335 

0.4695 

2.3 

0.0805 

0.0396 

-0.2885 

0.4302 

2.4 

0.0820 

-0.0081 

-0.2477 

0.3851 

2.5 

0.0792 

-0.0455 

-0.2116 

0.3364 
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p 

f(p) 

f'  (p) 

g(p) 

g'(p) 

2.6 

0.0732 

-0.0726 

-0.1804 

0.2869 

2.7 

0.0650 

-0.0898 

-0.1542 

0.2382 

2.8 

0.0556 

-0.0978 

-0.1327 

0.1926 

2.9 

0.0457 

-0.0982 

-0.1155 

0.1516 

3.0 

0.0361 

-0.0927 

-0.1022 

0.1161 

3.1 

0.0273 

-0.0828 

-0.0921 

0.0865 

3.2 

0.0197 

-0.0701 

-0.0846 

0.0634 

3.3 

0.0133 

-0.0566 

-0.0792 

0.0462 

3.4 

0.0084 

-0.0431 

-0.0752 

0.0341 

3.5 

0.0047 

-0.0307 

-0.0722 

0.0266 

3.6 

0.0021 

-0.0203 

-0.0698 

0.0224 

3.7 

0.0006 

-0.0118 

-0.0677 

0.0205 

3.8 

-0.0003 

-0.0056 

-0.0656 

0.0206 

3.9 

-0.0006 

-0.0013 

-0.0635 

0.0212 

4.0 

-0.0006 

0.0014 

-0.0614 

0.0223 

4.1 

-0.0004 

0.0024 

-0.0591 

0.0232 

4.2 

-0.0001 

0.0028 

-0.0567 

0.0236 

4.3 

0.0001 

0.0024 

-0.0544 

0.0237 

4.4 

0.0003 

0.0018 

-0.0520 

0.0232 

4.5 

0.0005 

0.0011 

-0.0497 

0.0224 

4.6 

0.0006 

0.0004 

-0.0476 

0.0213 

4.7 

0.0006 

-0.0001 

-0.0454 

0.0201 

4.8 

0.0005 

-0.0005 

-0.0435 

0.0188 

4.9 

0.0003 

-0.0008 

-0.0417 

0.0177 

5.0 

0.0004 

-0.0008 

-0.0400 

0.0166 

